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Abstract 

We present a detailed study of the electric dipole moment of the neutron induced by a vacuum 
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precision. With the current experimental bound this translates into a limit on the theta parameter 
of 1^1 < 3 X 10"^*^. We compare this result with the long-standing estimates obtained within chiral 
perturbation theory, and observe a numerical similarity, but also significant differences in the source 
of the dominant contribution. 
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I. INTRODUCTION 



Electromagnetic observables which are odd under T transformations are an important 
source of information about CP properties of the physics at and above the electroweak 
scale, complementary to that coming from K and B meson physics. In particular, there 
are now impressive experimental limits on the electric dipole moments (EDMs) of neutrons, 
heavy atoms, and molecules [H Q, 0, B 13 • The Kobayashi-Maskawa model, so successful 
in explaining the observed CP violation in K mesons, predicts EDMs to be several orders 
of magnitude smaller than the current experimental sensitivity. This presents a unique 
opportunity for limiting extra sources of CP-violation, and the constraints resulting from 
EDM data are generally very strong j^. 

In principle, EDMs can be used to probe the physics at a high energy scale by limiting 
the coefficients of operators Oi with dimension A; > 4 in the effective low energy Lagrangian. 
The effective Lagrangian for these operators has the form, 

£ ~ ^QM^-'^Of, (1) 

i 

where M is the mass scale at which these effective operators are induced and Cj their coeffi- 
cients which, in general, have logarithmic scale dependence. These operators are odd under 
CP transformations and their coefficients q are proportional to fundamental CP-violating 
phases of the underlying theory. 

Of these contributions, the electron EDM operator is the only example which may be 
constrained while avoiding the uncertainties necessarily associated with strong interactions. 
For the EDM of the neutron and the ^^^Hg atom, many more operators provide important 
contributions 0, 0, H, HO, IHl • Particular operators of interest include schematically, 

GG, ^qF(^lbq, j^qGa-fr,q, -^GGG, ^gPiggEzg, etc., (2) 

where F, G and q stand for electromagnetic, gluon and light quark fields, and Fi 2 denotes 
various contributing matrix structures l8|. 

In principle, the experiments 0, 0| impose strong constraints on the coefficients Cj. 
However, in practice, while the operators can be perturbatively evolved down to a scale 
of order 1 GeV, the ultimate connection between high energy parameters and low energy 
EDM observables necessarily involves non-perturbative physics. It is this final link which 
we wish to consider in this paper. The connection between different EDM observables and 
the coefficients q are especially important in supersymmetric theories where the number 
of relevant operators is much smaller than in a generic scenario and Cj can be explicitly 
calculated as a function of fundamental CP-violating SUSY phases (For a recent discussion 
in the context of the MSSM, see e.g. [ij). 

In the present paper, we focus on the first of these contributions, GG. This has a distinct 
status in that it has dimension=4 and thus receives contributions at tree-level from the 
fundamental QCD vacuum angle 6, the parameter labeling different super-selection sectors 
for the QCD vacuum. Experimental tests of CP symmetry suggest that 6 is small and, 
among different CP- violating observables, the EDM of the neutron {dn{0)) is the most 
sensitive to its value [l[ However, the calculation of dn{0) is a long standing problem 
[isi According to Ref. 0|, an estimate can be obtained within chiral perturbation 
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theory, relying on the numerical dominance of a one-loop diagram proportional to Inm^ 
near the chiral limit. The result can be conveniently expressed in the form 

and is seemingly justified near the chiral limit where the logarithmic term may dominate over 
other possible contributions parametrized in this formula by the constant c. However, these 
incalculable non-logarithmic contributions, can in principle be numerically more important 
than the logarithmic piece away from the chiral limit. In fact it is also worth noting that 
in the limit rriu, rrid 0, the logarithm is still finite, and stabilized, for example, by 
the electromagnetic mass difference between the proton and neutron. Consequently, one is 
unable to estimate the uncertainty of the prediction We note in passing that there 

is actually an additional 0(1) prefactor (1 — m^/m^) associated with the vanishing of the 
f/(l) anomaly at large N^.., which will play a role subsequently. The explicit derivation of 
this factor will be given in this paper. 

If the logarithm is cut off at the neutron mass, A ~ m„, and the non- logarithmic terms 
are ignored, one can derive a bound on the value of 9 using the current experimental results 
on the EDM of the neutron [113: ^ < 0(10-1°). Confronted with a naive expectation of 
^ ~ 1, the experimental evidence for a small if not zero value for 9 constitutes a serious fine 
tuning dilemma, usually referred to as the strong CP problem. As a consequence, one is 
usually led to introduce some mechanism via which the primary source of 9, the fundamental 
vacuum angle, is removed. However, even if this can be achieved, additional corrections are 
induced via the integration over heavy fields. Within this framework there are two main 
motivations for refining the calculation of dn{9). 

The first refers to theories where the axion mechanism is absent and the ^-parameter is 



zero at tree level as a result of exact P or CP symmetries jl5l . Il6l | . At a certain mass scale 
these symmetries are spontaneously broken and a nonzero 9 is induced through radiative 
corrections. At low energies, a radiatively induced theta term is the main source for the 
EDM of the neutron as other, higher dimensional, operators are negligibly small. As 9 itself 
can be reliably calculated when the model is specified, the main uncertainty in predicting 
the EDM comes from the calculation of dn{9). 

The second, and perhaps overriding, incentive to refine the calculation of dn{9) is due 
to efforts to limit CP-violating phases in supersymmetric theories in general, and in the 
Minimal Supersymmetric Standard Model (MSSM) in particular. Substantial CP-violating 
SUSY phases contribute significantly to 9 and therefore these models apparently require the 
existence of the axion mechanism. However, as mentioned above, this does not mean that 
the ^-parameter is identically zero. While removing 5 ~ 1, the axion vacuum will adjust 
itself to the minimum dictated by the presence of higher dimensional CP- violating operators 
which generate terms in the axionic potential linear in 9. This induced ^-parameter is then 
given by: 



9^nduced = where K^ = ^!. [ dxe''^^{0\T{^GG{x),Ocp{0m 



K \ ' ' 'Stt 



q=0 



where Ocp{0) can be any CP- violating operator with dim>4 composed from quark and 
gluon fields (as in Q), while 

K = zydxe^^^O\Ti^GG{x),^GGiom] ^ (4) 
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is the topological susceptibility correlator. In the case of the MSSM, the most important 
operators of this kind are colour electric dipole moments of light quarks qf^Gl^^^a^uInQ, 
and three-gluon CP-violating operators. The topological susceptibility correlator K was 
calculated in 17, 3| and the value of 9 generated by color EDMs can be found in a similar 



way 113] • Numerically, the contribution to the neutron EDM, arising from Oinduced is of the 
same order as direct contributions mediated by these operators and by the EDMs of quarks. 
Therefore, the complete calculation of dn as a function of the SUSY CP-violating phases 
must include a dn{0) contribution and a computation of this value, beyond the logarithmic 
estimate Q, is needed. 

In this paper, we present a detailed application of the QCD sum rule method ji^l to 
obtain an estimate for dn{9) beyond chiral perturbation theory. Within currently available 
analytic techniques, QCD sum rules seems the most promising approach to this problem as 
it has, in particular, been used successfully in the calculation of certain baryonic electro- 
magnetic form factors 21, i^l- Within the sum rule formalism, physical properties of the 



hadronic resonances are expressed via a combination of perturbative and nonperturbative 
contributions, the latter parametrized in terms of vacuum quark and gluon condensates. We 
note that previously QCD sum rules were used to estimate the neutron EDM induced by 
a CP-odd color electric dipole moment of quarks jl,!!^. Surprisingly, these results give dn 
~ 20 times smaller than the estimates based on the chiral loop approach The calcu- 
lation of dn{9) using QCD sum rules will certainly help to resolve this controversy. This 
question is of great numerical importance for the MSSM where color EDMs of quarks are 
large. 

The approach we shall use follows recent work i^j on the ^-induced p-meson EDM 



in reducing the operator product expansion to a set of vacuum condensates taken in an 
electromagnetic and topologically nontrivial background. Expansion to first order in 9 results 
in the appearance of matrix elements which can be calculated via the use of current algebra 
jiri Il8| . In this approach the 6'-dependence naturally arises with the correct quark mass 
dependence, and the relation to the U(l) problem becomes explicit as dn{9) vanishes when 
the mass of the U(l) "Goldstone boson" is set equal to the mass of pion. 

The initial results from this study were presented in 2^, and in this paper we shall 



present the details of this analysis. We also consider the relation of this result to the 
estimate obtained within chiral perturbation theory, and compare it with the outcome 
of an independent calculation of dn{9) reported recently in Ref. [26|. We begin in Section II 
by studying the phenomenological structure of the neutron correlator, and in particular 
addressing the issue of how to ensure chiral invariance of the result. In Section III we 
perform a tree level OPE analysis of the correlator to next-to-next-to leading order which 
corresponds to sensitivity at the level of 0{l/q^) terms. This requires an investigation of 
mixing with an additional set of currents CP-conjugate to the usual neutron interpolators. 
In Section IV, we combine the results of the previous two sections to construct a sum rule 
for dn{9) which we analyze numerically and extract the estimate, 

= 2.4 ± 1.0 X lO^^^^e ■ cm. (5) 

In section V we turn to chiral perturbation theory, and demonstrate explicitly the m^- 
dependence of the EDM and the CP-odd pion-nucleon coupling constant, indicating their 
explicit connection to the U{1) problem. We use this result to analyze and contrast the 
large N^. behavior of the chiral logarithm estimate and the QCD sum rule calculation of the 
EDM. We point out that, although unimportant numerically, the two results differ in the 
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large Nc limit, with the chiral estimate suppressed by a relative factor of 0{1/Nc). We then 
conclude in Section VI with some additional remarks. 



II. PHENOMENOLOGICAL PARAMETRIZATION AND CHIRAL INVARI- 
ANCE 

The starting point for the calculation is the correlator of currents rin{x) with quantum 
numbers of the neutron in a background with nonzero 6 and an electromagnetic field 

n(g2) = , J rfW«-(0|T{r/„(x)fJ„(0)}|0),,^, (6) 

where = —q"^, with q the current momentum. 

Before turning to the OPE analysis of this correlator, it is convenient to select an appro- 
priate Lorentz structure to consider and, in the present context, an important criterion will 
be invariance under chiral rotations. It is crucial to consider this issue when CP-symmetry 
is broken by a generic quark-gluon CP-violating source - the 6'-term in our case - as the 
coupling between the physical state (neutron) described by a spinor v and the current rjn 
then acquires an additional phase factor 

{0\vn\N) = XU^v, = e*"^^/2_ (7) 

The existence of this unphysical phase a is already apparent when one considers the sum 
rule for the neutron mass, which in the absence of CP-invariance can have an additional 
Dirac structure proportional to Z75. When we turn to electromagnetic form factors, this 
angle can mix electric {d) and magnetic (/i) dipole moment structures and complicate the 
extraction of d from the sum rule. 

To see how this will work, we recall that when considering 11 in the presence of some 
external field the phenomenological side of the sum rule may be parametrised by considering 
the form-factor Lagrangian which encodes the effective (in our case CP violating) vertices 
(see Fig. 1). We recall that after expanding to leading order in the background field, 11 is 
effectively a three-point correlator, and thus, although one can certainly write a two variable 
dispersion relation j21J], it lacks the powerful positivity constraints which follow from the 
analytic structure of the two-point correlator. Therefore, its more appropriate to instead 
explicitly parametrise 11 itself, rather than its discontinuity. The corresponding form-factor 
Lagrangian has the form C = fn.^il)^nS{q), where is the form factor, S{q) is the 
on-shell propagator for the neutron or one of its excited states, and (9„ is the operator 
corresponding to the induced vertex. 

Returning to the issue of chiral transformations, we can now consider the effect of such a 
mapping on the the double pole contribution on the phenomenological side of the sum rule. 
If we consider both electric and magnetic dipole moments, then the double pole term will 
have the form 

2(g^ -miy " 2iq^'^.j. U+^n){f^^a-dFa)U+m^). (8) 

in which we have introduced the dual field strength, F^^ = ^e^^apF'^^- Under a chiral 
rotation, we find that the numerator P transforms as follows 

UaPUa = mn{^, fJ'Fa — dFa} + m^{jiFa — (d + ajj,){Fa)) 

q'^ifiFa -{d- afi){Fa)) + 'iqf.quCr.yxi^F^u - {d - an)F^^), (9) 
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N • N' 

FIG. 1: Hadroiiic contributions to the current correlator in an external electromagnetic field. Possible 
excited states with the neutron quantum numbers are denoted generically by N' . 

where we have retained only the zeroth and first order terms of the expansion in a, and 
also neglected contributions proportional to ad. We see that only Lorentz structures with 
an odd number of 7-matrices are independent of a. In calculating dn, it is then clear that 
we should study the operator {Fa, as this is the unique choice with an unambiguous 
coefficient. 

The phenomenological side of the sum rule will therefore be parametrised in the form 

where, since we work outside the dispersion relation we may add polynomials in to ensure 
transversality in the chiral limit, and optimum behaviour for large q^. We then find that 
the function /(g^) takes the usual form, 

fi I \ ^ fij 1 1 ^ 

{q'^ — ml)"^ {q^ — rn^{iq^ — mf) ^ {q^ — rn'f){(f' — rn^Y 

where A is the coupling of the current to the neutron state ((Tj), (i„ is the neutron EDM, and 
fi and fij correspond respectively to transitions between the neutron and excited states, and 
between the excited states themselves. 

To suppress the contribution of excited states, we apply a Borel transform to 11, which 
we define, following |27|,i28|J, as 

f dY 

BU = lim3^„^oo,^/n=M2 ^^_-^y_ \ds) ^"^^^ 

where s = —q^. The continuum contributions in (jlip are then exponentially suppressed by a 
factor corresponding to the gap between and the next excited state, usually taken around 
(1.5GeV)^. However, while this suppression is quite large, previous studies of CP-even sum 
rules have found that the continuum contribution is not negligible (see e.g. jl^), and it is 
usually included for this reason. However, when studying correlators in background fields, as 
mentioned above, one is effectively dealing with a three-point function and one consequently 
loses positivity constraints for the contributions. Thus the couplings fij for example are not 
sign-definite (as would be the case for the two-point function). For this reason it seems 
inconsistent to parametrise the continuum in the normal way, and we have no alternative 
but to neglect it. In practice, we shall find that the sum rule we obtain is stable in any case. 

The coefficients of the single pole terms fi are also ambiguous in sign for the same reason, 
but these contributions are not exponentially suppressed by the Borel transform, and must 
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(a) 




(b) (c) 



FIG. 2: Various contributions to the CP-odd structure {Fa, 4}- (a) is the leading order contribution while 
(b) and (c) contribute at subleading order. 

be included for consistency. We then find that the phenomenological expression takes the 
form 

2 ' \{(f — — /' 

where the constant A parametrizes all the single pole contributions and, as we have ex- 
plained, is not sign definite. It is this expression that we shall contrast with the OPE 
calculation to be presented in the next section. 

III. CALCULATION OF THE WILSON OPE COEFFICIENTS 

We now turn to a tree-level calculation of 11 (jH)) within the framework of the operator 
product expansion. The neutron interpolating current Tjn is conveniently parametrised in 
the form. 

Tin = Jl+Pj2, (14) 

where the two contributions are given by 

ji = 2eabcid'^C'j5Ub)dc (15) 

j2 = 2tabc{dlCUb)^^dc. (16) 

The current j2 vanishes in the nonrelativistic limit, and lattice simulations have shown that 
ji indeed provides the dominant projection onto the neutron (see also ^3])- From 
we may define the coupling to the neutron state in the form (0|?7„|A^) = (Ai + /3A2)f , where 
As < Ai. 

Physical observables are independent of the parameter /3, and in principle any combina- 
tion of currents with neutron quantum numbers is appropriate. However, within the sum 
rules formalism, one has the imperative of suppressing the contribution of excited states and 



7 



higher dimensional operators in the OPE, and thus its convenient to choose (3 to this end. 
loffe has argued (sil. Is^. Issj ] that ~ —1 is an apparently optimal choice from a naly sis of 



the the mass sum rule (see also [29|). An argument based on minimal sensitivity [3J] leads 
instead to /3 ~ —0.2 using the same sum rule. We shall return to this issue again later, 
but for the moment it will be convenient to keep (3 arbitrary, and optimize once we have 
knowledge of the structure of the sum rule. 

In the absence of CP violation, ji and j2 form a basis for projection onto the neutron 
state. However, the presence of a CP violating source means that, in principle, it's also 
necessary to consider mixing with a CP-conjugate set of currents, which we shall write as 
follows, 

11 = 2eabc{d^Cub)dc (17) 

12 = '2eabc{dal5Cub)'y5dc- (18) 

However, we shall show subsequently that it's actually possible to choose a basis for 6 in 
which these two sets of currents do not mix, and the currents (ji, ^2) are a diagonalised 
combination for projection onto the neutron state. Thus for the time being we shall focus 
on T]n as the full current for calculation of H. 

We now proceed to study the OPE associated with (0). The relevant diagrams we need 
to consider are shown in Fig. 2 ((a), (b) and (c)). In parametrizing 6, we shall take a general 
initial condition in which a chiral rotation has been used to generate a 75-mass, so that 

C 5^ g^z75g/ + ea^G%G% + ■ ■ ■ , (19) 

in which we restrict to g/ = u, d, and so the reduced mass, which plays an important role in 
CP-odd observables, has the form, 

= 20 

m„ + vrtd 

The physical ^-parameter is of course = dq + Og-, but we shall keep the general form ()19|) 
and calculate the OPE as a function of both phases. The independence of the final answer 
from 9q — 9g will provide a nontrivial check on the consistency of our approach. We shall 
find that this requires the consideration of mixing with the additional currents (zi,i2), a 
point we shall come to shortly. 



A. Leading Order Contribution 

The leading order contribution is determined by the 1-loop diagram in Fig. 2(a). We work 
as usual with a constant background electromagnetic field, so that A^{x) = —^Ff^^{0)x'' , 
and for later reference we also use a fixed point gauge jl^ for the gluon potential, A^)f:"(x) = 
— ^G'^j^{0)t"'x'^ . At leading order it is not necessary to expand the quark propagator in the 
background field, but it is necessary to consider the short distance expansion of the quark 
wavefunction, 

q{x) = q{0) + XaD^q{0) + ■ ■ ■ , (21) 
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where Da = da — ieAa is the covariant derivative in the background field. When sandwiched 
between vacuum states, the second term contributes at first order in the quark mass via use 
of the equation of motion, pq = —imq. 

The vacuum structure is conveniently encoded in a generalized propagator which incor- 
porates the associated condensates. Projections onto particular vacuum condensates are 
chosen in order to obtain the Lorentz structure of interest (see e.g. 2^ for more details in 



the present context). The electromagnetic field dependence is determined in terms of the 
magnetic susceptibilities %, k and S,, introduced in j2l|: 

(0|ga^.g|0),. = xMmq\^) 
g{mGlX)Q.\^)F = sF^.(0|gg|0) (22) 
2^7(0|g75(GLt'^)g|0)F = ii,F,Mm\^). 

while the ^-dependence is either explicit in the case of 9q, or extracted via use of the 
anomalous Ward identity (see e.g. [l8|) in the case of Oq- This use of the anomalous Ward 
identity was discussed for example in [i^, and here we simply recall that the resulting 
expression for a generic structure mq{qTq)g^ has the form, 

mg{0\qrq\0)ea = ^m,^G(0|gr75g|0) + 0(mj), (23) 

where the ability to neglect the 0{im?) corrections follows since ^ m^r. The overall factor 



has the form (1 — m'^/m'i) jl8| which vanishes when t/(l)-symmetry is restored 



2J| . Making use of the anomaly, we see that the result then has formally the same form as 
arises from the 6q 75-mass contribution. This is once again a consequence of the anomaly, 
but an important point is that the sign of these contributions may differ, and thus we may 
obtain contributions having the unphysical form Oq — dq- The resolution of this puzzle will 
be described shortly. 

Defining iS{q) = {0\qa{x)qf^{0)\0) p^e, and ignoring a trivial 5-function over colour indices, 
the leading order propagator adapted to the CP-odd sector and appropriate for Fig. 2(a), 
then takes the form, 

-^^mFa^Xaiipl.U + ^(gg)(Fa(l + t9Gi,)U, (24) 

We shall henceforth follow [21] and assume that Xg = X^g stc, with fiavour independent 
parameters X: 

Substituting the generalized propagator into (jHI) according to the allowed contractions, 
performing the rather lengthy but straightforward algebraic manipulations, and Fourier 
transforming to momentum space, we find the result, 

A 2 

nLo(g') = -^(gg){Fa,^}[^(/? + l)2(4e,-e„) + 0~(l-/52)(2e,-ejlln^, (25) 
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where 9 = 9q + 9g, and 9 = 9g — 9q is an unphysical combination. The appearance of this 
unphysical combination is somewhat surprising and, as one might anticipate, is due to an 
additional source of mixing. The additional currents one needs to consider are precisely the 
CP conjugate set {ii,i2) introduced earlier on. 
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In order to illustrate the effect of this mixing we shall consider first the OPE for the 
two-point functions 

(0|T{jiJi}|0), (0|T{zi,zi}|0), (0|T{ji,zi}|0), (0|T{zi Ji}|0). (26) 

and the Lorentz structure proportional to ^. When both phases are set to zero, Oq = dq = 0, 
only the diagonal, ii and jj, two-point functions survive and the cross-terms in ()26|) are 
identically zero. When the 6'-phases are not zero, the ij correlators no longer vanish and an 
explicit calculation gives, 

I j d'xe^'^-^{Q\T{n;i,m = i{9G-eq)^^m,{qq) ^In^ (27) 
I f d'xe''^-mO\T{j,J,}\0) - {0\T{z,,H}m = T^m,{qq) ^In^. 

J ZTT^ — 

A straightforward re-diagonalization produces two "eigencurrents" , ji + ^ii and ii — ^ji, not 
mixed in the presence of 6g and dq. A similar procedure can be performed for the currents j2 
and i2- The mixing between ji and i2, and j2 and ii, is absent even for nonzero 9g, Og because 
these currents differ by an overall 75 matrix which effectively gives an anticommutator with 
q'. Therefore, the generalization of the neutron interpolating current (|T^ can be written in 
the following form, 

i9 

Vn = Jl + PJ2 + y (^1 + P^2) (28) 

Since the mixing between these two sets is explicitly proportional to Oc — Og, it is convenient 
to take 6q = as a useful choice of basis when working with ji and j2, where the mixing 
between ji 2 and ii^2, is simply absent. This situation resembles the problem in obtaining 
the "correct" quark mass behavior for the EDM of p, addressed in |24], in which one can 
simplify calculations by choosing m„ = rrid- Alternatively, one can use the generalized form 
for the current and observe that the presence of extra terms in ()28|) gives a contribution to 
the OPE cancelling identically the 6'-dependence of eq. (j^ . 

The outcome, either by including the mixing terms with {ii,i2), or by choosing the basis 
Gq = (^G) is the same, and we find 

A 2 

nLo(g') = 4} [^(/3 + 1)^(46, - e„)] In — ^, (29) 

explicitly proportional to the physical combination 9. 



B. Subleading Contributions 

Diagrams (b) and (c) in Fig. 2 require, in addition, the leading order expansion in the 
background gluon and electromagnetic fields, and S^lq is consequently more involved. We 
concentrate first on Fig. 2(b), which actually does not provide many non-zero contribu- 
tions. In addition to the leading order propagator presented in ()24|). the expansion in the 
background field leads to, 

Snlo = ^^^a/3(7/375)afe-^^^ln(-x2)(F(T(l-i^g75))„b 

+^{qq){i + ^OGi^)ab. (30) 
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With the field strength appearing exphcitly, it is not necessary to expand the quark wave- 
function, thus simphfying the calculation. The resulting expression for 11 actually receives no 
contributions from the first term of Snlo as it leads to a different Lorentz structure. How- 
ever, we see that since the propagator is logarithmic in the second term, the corresponding 
result behaves like ln(— which gives an infrared divergence in the Fourier transform 
to momentum space. These divergences were also observed in j2^, and we cut it off at a 
scale /i/ij. Strictly, the full OPE should be independent of such a cutoff, this dependence in 
the coefficient being cancelled by a similar dependence in the condensates. We shall ignore 
this subtlety here as these logarithmic terms will not enter our final sum rule. With this 
caveat, the resulting contribution to the {Fa, ^} structure in 11 has the form, 

m')^lo = ^(^^){^^' mP-lfe, (ln|^ - l) ^, (31) 

where we have checked that on calculating the mixing with (ii, ^2) the unphysical dependence 
on 6g — dq drops out and we are left with the expression presented here. 

Diagram (c) in Fig. 2 requires considerably more work, as there are several classes of 
contributions arising from it. Firstly, we need to expand the quark propagator in the back- 
ground gluon field, which gives rise to terms analogous to those for S^^q given above, since 
we make use of a fixed point gauge. Secondly, it is also possible to project such terms 
onto the vacuum in order to extract the leading dependence on F^^ using (j22I)- This re- 
quires a first order expansion of the quark wavefunction^ . The resulting propagator, ignoring 
contributions which lead to other Lorentz structures, may be written in the form 

where S"^^-^ is the expansion in the background gluon field. 



^& = -^,^G^p{l,l.)a, - ^ ln(-x^)(Ga(l - ^e,^,)U (33) 



while S^2) conveniently encodes the dependence on the condensates in the form. 



In calculating the corresponding contributions to 11, it is understood that one picks out the 
appropriate cross-terms in the correlator, in order to extract the leading order dependence 
on 9, and F^^. Infrared divergent logarithms also arise in this case from the final term in 
the propagator, while the other contributions actually contribute at next-to-next-to leading 
order, which in this case is 0{1/Q^). After a lengthy calculation, we obtain the following 
contributions to 11, 



i ((4/3^ - 4/3 + 2)e, + {3(3' + 2(3 + l)e„) 



(35) 



^ Note also that a second order expansion of the wavefunction, making use of the leading order 
propagator, would also provide contributions at subleading order, but these are combinatorially 
highly suppressed |2]| relative to the terms we are considering, and so will be ignored. 
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FIG. 3: Additional subleading contributions to the OPE. 



where we have checked that the unphysical 6g — Oq logarithmic terms are cancelled by mixing 
with (ii,i2), while we have simply used the "gauge" 9g = Oq to evaluate the 0(1/ Q"^) terms 
as these will turn out in fact to be numerically insignificant. 

Putting all the pieces together, we present the final result for the OPE structure arising 
from diagrams (a), (b), and (c) in Fig. 2, 



Is? 



X{(3 + l)\Aed - tu) In 



64772 

-4(/^-l)^e.(l + i(2. + 0) (l-J-l)^ 
-| ((4/32 - 4/5 + 2)erf + (3/32 + 2/5 + l)e„) ■ ■ ■ 



(36) 



This is our final result for the OPE, and we shall analyze the resulting sum rule obtained 
by equating this result with the phenomenological parametrization discussed in Section 2, 
in the next section. 

However, before turning to this analysis, we shall first make some comments regarding 
additional contributions that we have ignored. Some additional classes of diagrams are 
shown in Fig. 3. 

Naively, the diagram in Fig. 3(a) is loop suppressed relative to the contributions we 
have considered. This loop suppression is, however, fictitious as this diagram is pro- 
portional to the vacuum correlator as(47r)~^J (i^x(0|(G'(j'), (GG)|0) which is the same as 
J (i^x(0|(G'(j'), m^,gg|0). Nonetheless, we find in practice that all such contributions vanish. 

Diagrams of the form shown in Fig. 3(b) are more problematic because they involve 
correlators of the form 

~ / rf^xe^''^(0|T{(gg)2(x),m,5^g^z75g/}|0). (37) 

which are not calculable within our chiral approach. However, one suspects that these 
contributions of 0((gg)^), although suffering no loop-factor suppression, are small due in 
part to the small numerical size of (gg), but also due to combinatorial factors. We estimate 
these contributions via saturation with the physical t] meson. The result indeed turns out 
to be parametrically smaller than any term listed in Eq. ()36|1 . 

We shall conclude this section with some brief remarks on an independent calculation 
of dniff) using QCD sum rules, recently reported in [2^. The approach taken in this work 
is somewhat complementary to ours, as the authors use different Lorentz structures in @. 
These structures are not chirally invariant and have an admixture of the magnetic moment 
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a/i. Consequently, this requires simultaneous treatment of various two- and three-point 
correlation functions. The authors of introduce a nice chirally covariant notation to 



assist in keeping track of these terms. In their approach, one must then try to isolate the 
contributions to dn within sum rules that depend also on the neutron mass, /i, and the phase 
a. The subtlety here is that one must carefully keep track of all terms of 0{mg), despite the 
fact that fi for example is determined to leading order at 0(m°). As a result, the separation 
of c? ~ m6{qq) arises as a delicate cancelation between a combination of different terms, each 
of order 0(m°). The authors of [2^ simplify matters somewhat by setting all quark masses 
equal rriu = riid, and therefore don't observe the appearance of the reduced mass m* (plj) . 

_At first sight there is a significant problem with the use of Lorentz structures other than 
{Fa,c^, which is the appearance at leading order of certain incalculable condensates of the 
form 

(0|r{gz75g(a:),m* ^gyZ75g/}|0) and (0|T{ga^^g(x), ^ gyi75g/}|0). (38) 
/ / 

In our approach these terms arise manifestly at 0{mg), while for other channels these terms 
can arise at 0{mq) and apparently need to be dealt with. The first condensate in ()38p 
can be connected to an 0(m*) correction to the correlator J d'^x{0\T{qij5q{x), {GG){0)}\0) , 
for which no reliable means of extracting its value is known (for a recent discussion of the 
subleading mass dependence of this and related correlators see, e.g. [36]). These terms, as 
well as ^ and /t-proportional contributions, are ignored in ji^. However, this is consistent 
as one may show that, remarkably enough, in the final result for the ratio dn/lJ-n these 
corrections only appear at 0{mq{qq)'^); a subleading effect. This wasn't explicitly pointed 
out in [2^, but we believe it is an important (and indeed interesting) point, necessary for the 
consistency of their approach. Numerically, the results of jl^ are very close to the results 
we obtained in jl^, and shall discuss in the next section. Thus, the two approaches indeed 
appear quite complementary. 

As a final related comment, we note that the use of chirally non- invariant channels may 
be the origin of numerical problems in the QCD sum rule calculations of dn induced by color 
EDMs jsllisl. This, of course, is yet to be checked but our preliminary estimates show that 
the use of the chirally invariant channel yields the neutron EDM at a level comparable with 
phenomenological estimates 



IV. NUMERICAL ANALYSIS OF THE SUM RULE 

Turning now to the analysis of the sum rule ()36|). an inspection indicates that the standard 
choice of /3 = — 1, appropriate for CP-even sum rules, will not be the most appropriate here 
as it removes the leading order contribution. In general there are two motivated criteria for 
fixing the mixing parameter /5: (1) at a local extremum 1341 : or (2) to minimize the effects 
of the continuum and higher dimensional operators jiol l3lL l3^. l33|. We find in this case 
that extremizing in the parameter (3 also leads to the unappealing cancelation of the leading 
order contribution. Thus the most natural procedure appears to be to choose (3 in order 
to cancel the subleading infrared logarithm which is ambiguous as a result of the required 
infrared cutoff. This procedure actually mimics the effect of the choice (3 = —1 in the sum 
rule for the nucleon mass. We therefore take (3=1, and it is this choice that we shall now 
contrast with the phenomenological side of the sum rule. It is important to note, however. 
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that use of the "lattice" current with (3 = will also produce a numerically similar result. 
We shall make further comments on this issue at the end of the section. 

On the phenomenological side of the sum rule we have ()13|) in which the coupling is now 
interpreted as A = Ai + PX2. After a Borel transform of ipHj) and (fT^. and using /? = 1 as 
discussed above, we obtain the sum rule 



X rrindn + AM' 



327r2 
4x(4erf - e„) 



2M2 



e(4ed + Sen) 



(39) 



The coupling A present in fl39p may be obtained from the well known sum rules for the 
tensor structures 1 and ^ in the CP even sector (see e.g. j2£| for a recent review). 

To aid in the presentation of the results it is convenient to define an additional function 



z/(M2 



(40) 



which is then determined by the right hand side of (jHOI)- Inspection of (jHEI) suggests that it is 
not appropriate in this case to try and remove the unknown parameter A, via differentiation 
for example. Instead we shall make the conventional assumption that the left hand side is 
a linear function of (i.e. A is independent of M^), and construct two sum rules whose 
behaviour will allow us to estimate the slope of this line, and thus the parameter A. In fact 
this approach will lead to a result which for consistency will require A ~ 0, and thus allow 
an extraction of the EDM parameter dn- We now construct these two sum rules as follows: 

• (a) Firstly, we extract a numerical value for A via a direct analysis of the CP even 
sum rules. This analysis has been discussed before and will not be reproduced here 
(see e.g. jl^). One obtains^ 

(27r)^Aa ~ 1.05 ±0.1, (41) 
which leads to a sum rule of the form. 



{qq)e 



4x(4ed - eu) - ^^(4ed + 8e„ 



(42) 



(b) As an alternative, we extract A explicitly as a function of (3 from the CP-even sum 
rule for ^ [2_9] which we reproduce here for completeness. 



{2'KYX'e 



5 + 2(3 + 5(3^' 



64 



i/M2 



5 + 2(3 + 5(3"' 



256 



2M4 

-S2/M' 



M2 



+ 1 



(43) 



^ The particular choice of (5 in the range [—1, 1] is not very important here as the numerical value 
for A is not highly sensitive to this choice. 
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FIG. 4: The neutron EDM function i^(M^(GeV^)) is plotted according to the sum rules (a) and (b). The 
dashed line shows the contribution from the leading order term only. 

where Si and S2 parametrise the continuum thresholds, while h = (27r)^(^G'^) ~ 
0.47± 0.2GeV'^, and we have neglected higher dimensional contributions, and leading- 
log anomalous dimension factors. The reason for the omission of the latter is that 
they provide a negligible contribution when is small, as will be the case here, and 
furthermore at this scale one has good reason to distrust the leading log approximation. 
Therefore, the (estimated) effect of such factors will be combined into the overall error 
estimate. 

Solving for \h, we obtain a new CP-odd sum rule by substituting the result into 
dSH) and setting /3 = 1, 



''^ ' 647r2m„^^^^ 4M6ci(M2,si) + 6M2c2(M2,S2) ■ ^ ^ 

In this expression ci and C2 are the continuum parametrizations introduced in ()43j] . 
Throughout we shall assume si = S2- 

We shall now proceed to analyze the sum rules numerically. Note that our assumption 
that z/(M2) is linear in also requires that A is constant in an appropriate range of the 
Borel parameter. This point can be checked explicitly in case (b) above. 

The two sum rules described above for z/^ and Uf, are plotted in Fig. 4, where the effect 
of the higher dimensional terms in ()39p proportional to ^ is also displayed. //(M^) is to be 
interpreted as a tangent to the curves in Fig. 4. For numerical calculation we make use of 
the following parameter values: For the quark condensate, we take 

(0|gg|0) = -(0.225 GeV)^ (45) 

while for the condensate susceptibilities, we have the values 



X = -5.7 ±0.6 GeV"^ j37|, (46) 
^ = -0.74 ±0.2 Q. (47) 

Note that which enters at 0(1/M2), since it is dimensionful, is numerically significantly 
larger than ^. In extracting A in case (b) we also set a relatively large continuum threshold 
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So = (2GeV)^ for consistency with the CP-odd sum rule in which this continuum is ignored 
for reasons discussed earher. 

One observes that both sum rules have extrema consistent to ~ 10%, suggesting that our 
procedure for fixing the parameter j3 is appropriate. Furthermore, the differing behaviour 
away from the extrema implies that for consistency we must assume A to be small. One 
then finds dn ~ v{M'^ ~ 0.5GeV^). It is also interesting that the effective scale is around 
M ~ 0.7GeV which is well below m„, and should be cause for concern regarding the con- 
vergence of the OPE. Nonetheless, one sees that the corrections associated with the leading 
higher dimensional operators are still quite small. This low scale is also the reason we have 
ignored leading-log estimates for the anomalous dimensions as noted earlier in the context 
of extracting the coupling A, as their status is unclear at this scale. A naive application 
leads to a small correction that we shall subsume into our error estimate. 

Extracting a numerical estimate for dn from Fig. 2, and determining an approximate 
error arising from: (1) analysis of the sum rule; (2) the error in x; and (3) an estimate of 
ibO(20%) for higher dimensional operators and anomalous dimension factors; we find the 
result^ 

dn = (1.0 ± 0.4)e^— . = (3.6 ± 1.4) ■ lO'^eg , , (48) 

^ ' (500MeV)2 ^ ' (100MeV)5 (m„ + md)2' ^ ' 

for the neutron EDM, for which the dominant contribution naturally arises from x- Com- 
parison with the result of Ref. f]3] indicates rather good agreement in magnitude, due 
essentially to the low effective mass scale M ~ 700MeV. We also obtain dn of the same sign 
if one assumes no significant corrections to the logarithm in ^. The relation between 
the calculation presented here and the chiral logarithm estimate will be discussed in more 
detail in the next section. 

However, before turning to this comparison, we would like to comment further on the 
specific choice of the current, i.e. the choice of (3. As mentioned earlier, a conventional 
choice for CP-even sum rules is that advocated by loffe [33!]: (3 = —1. However, even in 
these channels, standard minimal sensitivity arguments, (as used for example in the context 
of renormalization scheme dependence jl^), lead instead to a choice of /3 ~ —0.2 jj^ ]. 
Furthermore, we have found that an apparently optimal choice for the chirally invariant 
CP-odd channel is instead [3 = 1. 

This variation in (3 might be interpreted as a sign of inherent uncertainties due to the 
effect of excited states or higher order condensates. However, we would like to point out that 
this need not be the case if these differences are observed for different physical observables, 
as is the case here. In particular, while one expects that the optimization of j3 with respect 
to minimizing contamination from excited states may be universal, loffe [i^ has emphasized 
that this aim conflicts with the need to minimize the uncertainties due to higher dimensional 
operators. The latter point is quite "observable dependent" , and thus the best compromise 
choice for (3 may differ for different observables. In particular, as pointed out in |23|, the 
points /? ~ — 1 and -|-1 are not distinguished purely on the basis of removing contamination 
by excited states. 



^ (04/2005) v4: This updated expression corrects an overall factor of two error in previous versions. 
^ The sign of dn{0) in differs simply because of an overall minus sign in the initial definition of 
the 0-term. 
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FIG. 5: Leading contribution as ^ to (i„ in the chiral loop approach There is also a second 

diagram with the pion-nucleon vertices interchanged. 

While we find such arguments for the overall consistency of the current with j3 = 1 
compelling enough, we have also exphcitly tested the /3-dependence of our result. To do 
this, it is necessary to consider the effect of the logarithmic terms in ()36|) . Cutting the logs 
at fiiR ~ ^QCDi we then find that the numerical value for does not vary that dramatically 
with (3. Indeed for /3 ~ 0, we find a result which is numerically within ~ 30% of the result 
obtained at /3 = 1. 



V. RELATION WITH ESTIMATES FROM CHIRAL PERTURBATION THEORY 

In this section, we shall consider the relation between the value for dn we have extracted 
using QCD sum rules, and the well known chiral loop calculation jT3]. Within the sum rules 
approach, use of the chiral anomaly made it quite explicit that the result would necessarily 
vanish in the limit m^r. To compare this with the chiral result, we first need to observe 

how the latter also satisfies this necessary consistency condition. As far as we are aware, 
this result has not appeared in the literature and we present a derivation of the expected 
prefactor (1 — m^/m^) in the next subsection. 

Prior to addressing this issue, we first recall the main points of the chiral approach 
used in . A chiral rotation is used to recast the 6'-term into a CP violating quark mass, 
5C.CP = iOm^: J2f Qfl^Qf- This induces an effective CP-odd interaction of pions with nucleons 
in terms of a small correction to the standard 75-interaction, 

jCnNN = TT''NT''{ig^NNl5 + 9nNN)N. (49) 

The CP-odd pion-nucleon coupling ^^jv^v induced by 9 may be obtained by the PCAC 
reduction of the pion 

g^NNNr'^N = -_(iV| [<5£cp, ™] |iV) = —^{N\qT''q\N). (50) 

^ Jtt Jtt 

The remaining matrix element can then be expressed in terms of the mass splittings in the 
baryon octet using SU (3) flavour symmetry. Note that the expression (jMIj) is valid only for 
soft pions with momenta smaller than the characteristic hadronic scale of 1 GeV. 

The crucial observation made in 3| is that the electric dipole moment of the neutron. 



induced by the source SCqp, can be estimated using the singular behaviour of the chiral 
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FIG. 6: Different contributions to the g-„NN coupling. The leading diagram (a) is identically cancelled by 
(b) and (c) in the limit — s- TOtt. 



loop shown in Fig. 5. The loop integral has a logarithmic infrared divergence as m,r 
and one then obtains the estimate [14] 

where the cutoff represents the scale at which chiral perturbation theory breaks down. It is 
reasonable to assume that A ~ mp or m„. After substituting the numerical expression for 
Q-kNN we arrive at the result reproduced in Q. 



A. C/A(l)-properties of the EDM in chiral loop calculations 

The independence of the original chiral loop result (j3T|) from m„ has been the source of 
some controversy, and continued work, in the literature (see e.g. j39|). Before we compare 
our results with that of Ref [l^l, we wish to point out the particular mechanism, ignored in 
the discussion above, which sets the EDM to zero when Ua{^) symmetry is restored. This 
cancelation should, of course, be very similar to the vanishing of the CP-odd amplitude for 
T] ^ TTTT, demonstrated in |l8|. 

Within the chiral loop approach reviewed above, the EDM depends explicitly on the 
CP-odd pion-nucleon coupling Qt^nn induced by 9. We shall now prove that this coupling 
is proportional to the factor (1 — m'^/m'^) for the case of two flavors, taking m„ = = ruq 
for simplicity. 

Recall that in the calculation of 3|; the 6'-term is written in the form of a singlet com- 



bination of quark 75-mass terms. Since the ?7-meson has the same quantum numbers, this 
combination can produce rj from vacuum with an amplitude proportional to f~^mq{0\qq\0). 
Thus, the calculation of (jnNN should account for the additional contributions related to r] be- 
ing produced from the vacuum and then reabsorbed by the nucleon, subsequently producing 
the soft pion (see Fig 6). 

In the absence of any physical effect from the anomaly, the mass of 77 should vanish in 
the chiral limit. In this case the fiavor-singlet field rj can be treated exactly as a pion field. 
Therefore, the amplitude for the low energy scattering tcN rjN can be related to the 
nucleon sigma term in the triplet channel, 

M^a^^^^ ~ mg{N\qr^q\N), (52) 
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which is exactly the structure appearing in the calculation of the CP-odd pion-nucleon 
coupling constant g-^NN, eq. (|MH). 

Summing different contributions we obtain 

,^,,Nr'N = p{NWr',\N) ( 1 - '?^) = ^(iV|gr»,|7V> (63) 

with the anticipated dependence on m^. 

Now the connection to the Ua^^) problem and the vanishing of Qt^nn in the limit 
become explicit. In real life, due to the physical effect of the anomaly, m,, remains finite in 
the chiral limit, whereas m^r — >■ 0, and the correction from diagrams (b) and (c) is negligible 
as it is second order in the light quark masses. Thus numerical extraction of bounds on 6 
using this result, as in are essentially unaffected by this correction. Nevertheless we 
find this exercise rather instructive, noting that the coupling gT^NNi^) may also be used for 
the extraction of the limit on 6 from the mercury EDM experiment (3|]. 

As a small digression, it is also worth pointing out that in where the issue of Ua{1) 
restoration was also addressed, it was suggested that contributions to dn{0) can be obtained 
within the chiral approach which are directly proportional to the anomalous magnetic mo- 
ment of neutron. This seems highly unlikely since fi and d are generated at quite different 
energy scales. The one-loop diagram contributing to d receives contributions from momen- 
tum scales ml <C <C IGeV^, whereas similar diagrams for fi diverge quadratically and 
thus are saturated by momenta ~ 1 GeV where the chiral description breaks down. 

B. Relation With the EDM Prom Chiral Perturbation Theory 

Given the numerical similarity between the result we have obtained for dn{0) and the 
estimates based on the dominance of the chiral logarithm, Inm^^, it is natural to ask whether 
or not this logarithm is hidden somewhere in the OPE analysis. 

Such a suggestion is necessarily speculative, as the calculations are performed using dif- 
ferent dynamical degrees of freedom, and are in principle valid at quite different external 
momentum scales. Nonetheless, at first sight, one may be tempted to identify the chiral log 
term explicitly with the subleading infrared log-terms obtained in (j36|) . The full momen- 
tum dependence of this term is, however, lii{Q^ / fj,'^) / Q"^ , and thus the singular behaviour is 
instead determined by Nevertheless, one might suggest that this naive extrapolation 

to the chiral regime is not appropriate, and the power-like term gets softened, while the log- 
arithm remains. The most obvious log-term of this kind arises from diagram (b) in Fig. 2. 
Assuming the top two quark lines are soft, as appropriate for the condensate, we can inter- 
pret this as a soft pion, and by adding two spectator quarks, we get something analogous 
to the chiral loop diagram contributing in the approach of [1J|. Here, CP-violation arises 
from a particular 4-fermion vertex proportional to 6. 

While such heuristic relations are sometimes possible when the momentum dependence is 
easily equated on both sides, we would now like to argue that, at least at large Nc, the chiral 
logarithm is actually subleading with respect to the leading order term in the OPE. While 
this need not be important numerically for Nc = 3, it provides a convenient parametric 
distinction between the results obtained using sum rules and chiral perturbation theory. 

At large Nc, the Ua{^) restoration factor for the chiral expression derived in the previous 
subsection becomes important. In particular, as is clear from the discussion above, the result 
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of calculation, either with the use of chiral techniques or QCD sum rules, is proportional 
to (1 — m^/m^), and thus the EDM has a natural parametric dependence on 9 of the form 

However, the main distinction between the sum rules and chiral loop calculations is that 
the quark diagrams used within the sum rules calculation are effectively tree-level. Thus 
one expects the effect to be parametrically enhanced relative to the chiral loop. Indeed, 
this follows from an additional suppression factor of 1//^, which is 0{1/N^, and which is 
absent in the OPE calculation. The corresponding factor has the form {qq)/mn, which is 
0{Nc/Nc ~ 1). Summarizing, we obtain the following behaviour at large in these two 
cases, 

QCD-SR oc P^^^0{9/N^) 
Chiral Loop d„ oc -dd\N) ^ ^^^/^.^ ^^^^ 

where we have only kept track of the relative A^c^dependent factors^ It is interesting to 
note that the isovector matrix element {N\uu — dd\N) contributing to the chiral result does 
not grow with Nc-, whereas both {N\uu\N) and {N\dd\N) are proportional to N^. This is 
because the neutron in the large A^^^ limit has {N^. — l)/2 w-quarks and [Nf. + l)/2 d-quarks. 

In summary, while there appears to be at least a qualitative mechanism for mapping 
some of the OPE diagrams to those which would produce a chiral logarithm, the behaviour 
at large is of the form (c + \n.m.,^/Nc). Of course, we reiterate that physically there need 
be no suppression for Nc = ?>. 



VI. DISCUSSION 

The calculation of dn{0) via QCD sum rules has produced a numerical result very close 
to the estimates obtained using different techniques. Indeed the chiral loop estimate, QCD 
sum rule calculations, and even the naive quark model [l^ , all agree and predict essentially 
the same EDM to within a factor of 2-3. It is interesting to note that the power counting 
procedure in the chiral theory, combined with the quark model for nucleons known as naive 



dimensional analysis |41[, also produces a similar estimate. 

The situation is apparently very different in the case of dimension 5 and 6 CP-odd 
operators, especially for color EDMs. Different methods produce results varying within 
more than one order of magnitude. Moreover, none of the existing calculations is capable 
of answering the question of which combination of color EDMs of m, d and s quarks in fact 
enters into the single observable dn- We believe that the approach developed in this paper 
can be used for the calculation of the neutron EDM induced by these dimension 5 and 6 
CP-odd operators, and such a calculation will clarify this issue. A technical outcome of the 
work presented here is that the analysis can be done in the chirally invariant channel, where 
the OPE is directly proportional to the electric dipole moment. 



^ An additional Nc enhancement, due to the scaling of qa ^ Nc has been discussed in relation to Skyrme 
model calculations (see e.g. This effect should be generic to both approaches, and we have ignored 

such overall contributions. 
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There is an additional incentive for using QCD sum rules for the calculation of dn, as 
opposed to any other method. In the case of dim=5 operators, QCD sum rules would 
normally produce in the form of a linear combination of terms like dq{qq), where dq is an 
EDM or a color EDM of a light quark. In many models, including the MSSM and variants 
thereof, dg ~ nig which would allow the removal of much of the uncertainty related to the 
imprecise knowledge of light quark masses. Indeed, when multiplied by the value of the 
quark condensate, any linear combination of dg can be expressed in terms of Tn^f^ times a 
function depending only on the ratio of light quark masses. The latter is known to much 
better accuracy than the quark masses themselves and does not depend on the normalization 
scale. 

In conclusion, we have presented a QCD sum rules calculation of the 6'-induced neutron 
EDM. This result is explicitly tied to a set of vacuum correlators which are non-vanishing 
only in the absence of a U(l) "Goldstone boson". The use of QCD sum rules in the chirally 
invariant channel allowed us to unambiguously extract dn{0), and independence of the answer 
from any particular representation of the theta term (jl9|) was checked explicitly. 

Combining our result with the recently improved experimental bound on dn [2] we derive 
the limit on theta: 

1^1 < 3 X 10"^°, (55) 

which is quite close to previous bounds, and actually somewhat less constraining as our result 
for dn{0) is slightly lower than the corresponding result obtained within chiral perturbation 
theory with the cutoff at rup. Numerically, to a large extent, our result is proportional to x, 
the electromagnetic susceptibility of the vacuum. 
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